PASJ: Publ. Astron. Soc. Japan ??, 1-??, 
© 2008. Astronomical Society of Japan. 



A Measurement of the Quadrupole Power Spectrum 
in the Clustering of the 2dF QSO Survey 



Kazuhiro Yamamoto, Masashi Nakamichi, Akinari Kamino 
Graduate School of Science, Hiroshima University, Higashi-Hiroshima, 735-8526, Japan 

Bruce A. Bassett 
Institute of Cosmology and Gravitation, University of Portsmouth, England and 
SAAO, Observatory, Cape Town, South Africa 
Hiroaki Nishioka 

Institute of Astronomy and Astrophysics, Academia Sinica, Taipei 106, Taiwan, R.O.C. 



We report a measurement of the quadrupole power spectrum in the two degree field (2dF) QSO rcdshift 
(2QZ) survey. The analysis uses an algorithm parallel to that for the estimation of the standard monopole 
power spectrum without first requiring computation of the correlation function or the anisotropic power 
spectrum. The error on the quadrupole spectrum is rather large but the best fit value of the bias parameter 
from the quadrupole spectrum is consistent with that from previous investigations of the 2dF data. 

Key words: methods: numerical - quasars: general - cosmological parameters - large-scale structure 
of Universe 



1. Introduction 

Since large cosmological surveys are carried out in redshift space, the peculiar velocity of sources distorts the apparent 
spatial distribution of cosmological objects. These redshift-space distortions are one of the important effects in large 
redshift surveys (e.g., Suto et al. 2000a, Magira, Jing, Suto 2000). Redshift distortions can be classified into the linear 
distortion arising in the linear theory of density perturbation and the Finger of God effect arising in the nonlinear 
regime. It has been also pointed out that the the geometrical effect from the expansion of the universe causes another 
apparent distortion in the distribution of cosmological objects, which is referred to as the cosmological rcdshift-spacc 
distortion or the geometric distortion (Alcock, Paczynski 1979, Ballinger, Peacock, Heavens 1996, Matsubara, Suto 
1996). Thus the detection of the redshift-space distortion provides us with unique information about the peculiar 
velocity and expansion history of the universe. 

In the Kilo-Aperture Optical Spectrograph (KAOS/WFMOS, see http://www.noao.edu/kaos/), one of the future 
large survey projects, it will be possible to confront in great detail the theoretical redshift-space distortions including 
the geometric distortion with observation (e.g., Yamamoto, Nishioka, Bassett 2005). Such surveys will provide signif- 
icant dark energy constraints (Matsubara, Szalay 2003, Seo, Eisenstein 2003, Lindcr 2003, Blake, Glazcbrook 2003, 
Yamamoto 2003;2004, Amcndola ct al.2005) using the baryon acoustic signatures, which have recently been detected 
clearly in the SDSS and 2df data sets. (Eisenstein et al. 2005, Cole et al. 2005, Yahata et al. 2005). 

In general it is known that the redshift-space power spectrum on a constant hypcrsurface of the redshift z can be 
expanded as (Taylor, Hamilton 1996) 



where Ci(fx) is the Legendre polynomials, /i(= cos (9) is the directional cosine between the line of sight direction and 
k (See Figure 1). 1 The monopole Po(k,z) represents the angular averaged power spectrum, what we mean by 
the power spectrum usually, and the quadrupole P2(k,z) represents the leading anisotropies in the power spectrum 
due to the redshift-space distortions. The odd moments vanish by symmetry. The quadrupole spectrum provides us 
with new information over what is available in the monopole power spectrum. In practice, the quadrupole spectrum 
reflects the peculiar and random velocity of samples (Kaiser 1987, Cole, Fisher, Weinberg 1994, Hamilton 1996). 
It has been shown that the quadrupole spectrum will be especially useful in breaking the degeneracy between the 
bias and the dark energy parameters when measurements of the power spectrum are good enough in future surveys 
(Yamamoto et al. 2005). Thus accurate measurement of the quadrupole spectrum can be important in characterizing 

1 Note that our definition of the multipole spectrum P; is different from the conventional definition by the factor 22 + 1. 
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the redshift-space distortions quantitatively. 

Pioneering work on the measurement of the quadrupole power spectrum was carried out by Cole, Fisher and 
Weinberg (1994) and Hamilton (1996, 1997) using galaxy rcdshift survey catalogs. Cole et al. (1994) presented a 
systematic method to estimate the quadrupole power spectrum through the anisotropic power spectrum. Using the 
method, Hatton and Cole (1999) estimated the /3-factor from the quadrupole power spectrum in the two degree field 
(2dF) galaxy survey. In the work by Hamilton (1996,1997) the quadrupole power spectrum was obtained by the 
transformation of the correlation function. In the present work, however, we consider a different method to estimate 
the higher multipole moments of the power spectrum. Our method is parallel to the one widely used to measure the 
monopole power spectrum (Feldman, Kaiser, Peacock 1994, Yamamoto 2003) and allows us to obtain the multipoles 
of the redshift-space power spectrum without evaluating the correlation function or the anisotropic power spectrum. 
The difference between our method and that developed by Cole ct al.(1994) is discussed in the next section. 

This paper is organized as follows: In section 2, we explain the algorithm for estimating the multipole moments of 
the power spectrum. In section 3, we apply the method to the 2dF QSO sample and obtain the quadrupole power 
spectrum. Using the result we briefly discuss the constraint on the bias parameter. The last section is devoted to 
summary and discussions. Throughout this paper we use units in which the speed of light is unity, c = 1. 



2. Method 



Here we explain the optimal weighting scheme for estimating the multipole moments of the anisotropic (i.e. redshift- 
space) power spectrum. Our method developed here is unique and different from those developed previously. Our 
method does not require computation of the correlation function; c.f. Hamilton (1996). A similar method was 
developed by Cole et al. (1994), which also does not require the computation of the correlation function. However 
it requires estimating the anisotropic power spectrum. The difference between our method and their method will be 
discussed in detail below, however, the essence is as follows. 

In the Cole et al. method, it is necessary to divide the sample into subsamples and introduce one line of sight 
direction for each subsample in order to explicitly define /z. Thus in the Cole et al method, the line of sight directions 
for all objects are regarded as the same in each divided subsample. In our method, however, we need neither divide 
the sample nor explicitly define fx, instead, we introduce fi for each pair of objects of the sample (see Figure 1). 
Furthermore, the multipole moments can be computed fully in parallel with the conventional monopole estimation 
(i.e. the 'usual power spectrum'). 

Throughout this paper, we use s to denote the three dimensional coordinates in redshift space. We also use n g (s) 
and n(s) to denote real catalog and the expected mean number density. A fluctuation field is then defined via 

F(s) = n g {s)-an s (s), (2) 

where n g (s) = ^ 5(s — s^) with being the location of the ith object of real catalog; similarly n s (s) — J^j $( s ~ s j) 
is the density of a synthetic catalog. The synthetic catalog is a set of random points without correlation, which can 
be constructed through random process by mimicking the selection function of the real catalog. Here we assume that 
the synthetic catalog has a mean number density 1/a times that of the real catalog. The synthetic catalog provides 
an estimate for the mean number density in the absence of clustering, then it is useful to provide an estimator for the 
fluctuation field F(s). From the definition of the random field n s (s), we assume (Feldman ct al. 1994), 

(n g (s 1 )n g (s 2 )) = n(si)n(s 2 )(l + £(si,s 2 ))+n(si)<5(si - s 2 ), (3) 

(n s (si)n s (s 2 )) = a^ 2 n(s 1 )n(s 2 ) + a~ 1 n(s 1 )6(s 1 - s 2 ), (4) 

(n g (si)n s (s 2 )) = a~ 1 n(si)fi(s 2 ), (5) 

where £(si,s 2 ) denotes the two-point correlation function. Using these relations, we have 

(F( Sl )F(s 2 )) = n(si)n(s 2 )£(si,s 2 ) + (1 + a)n( Sl )S( Sl - s 2 ). (6) 

The estimator for the multipole moments of the power spectrum may be defined as follows (see Yamamoto, Nishioka, 
Taruya 2000), 

R f,. = /ds 1 /ds 2 y>(s 1 ,k)V;(s 2 ,k)F(s 1 )F(s 2 )e' k ^- s ^A(§/ t -k) 
l[ ' Jdsn 2 (s)^(s,k) 2 ' 1 ' 

where Sh = (si + s 2 )/2 and the 'hat' means the unit vector, and ij}{ s M) is a weight factor. As mentioned before, £i(fi) 
is the Lcgendre polynomials. (See also Figure 1 for the definition of variables). We can choose the weight factor so as 
to minimize the variance of the spectrum in order to optimize the performance of the estimator according to the local 
density field (see below). 

Here we adopt the approximate formula 



No. ?] 3 

/rile 
— P( k} \ Sh \)e-^-^, (8) 

where -P(k, |s^|) is the power spectrum defined on a constant-time hypersurface. This expression is not strictly correct 
when the rcdshift Si and S2 are very different. Thus (8) is justified within the distant observer approximation, 
|si — 82 1 -C |si|, |s2 1 , i.e., in case of short separation between si and S2. 

Using the above relations and under the distant observer approximation, the ensemble average of Ri(k) is 

/rfsn(s)V(s,k) 2 P(k,|s|)A(s-k) 
(i?i(k)> " /dsn 2 (s)V>(s,k) 2 + 5i(k) ' (9) 

where 

(l + «)/rfsn(s)^s,k) 2 £ ; (s-k) 

/dsn 2 (s)<A(s,k) 2 ■ (iUj 

The estimator of the anisotropic power spectrum is obtained by subtracting the shotnoise Sj(k), yielding 

V l (k) = R l (k)-S l (k). (11) 
Next we focus on the covariance of Vi(k), which is given by 

(ATMk)AW)) = m(k)- (vMWiM- (W))}) 

= (V l (k)V l (k'))-(V l (k))(V l (k')) 

= (R l (k)R l (k>))-(R l (k))(R l (k>)). (12) 

Assuming Gaussian statistics some calculation gives: 

(F(s 1 )F(s 2 )F(s 3 ) J F 1 (s4)) = ( J F 1 (s 1 )F(s 2 ))(F(s 3 ) J F 1 (s4)) 
+ (P(s 1 )F(s 3 ))(F(s 2 )P(s 4 )) 

+ (F( Sl )F(s 4 ))(F(s 3 )F(s 2 )), (13) 

we have 

/ l,'\2l-l 



(A^(k)A^(k')) = [ J dsn(s)V(s,k) 2 ]- 1 [ J ds'n(s') 2 ^(s',k') 



2 4 

x2f[[J ds^^flij ds^(s J ,k')}(F(s 1 )F(s 3 ))(F(s 2 )F(s 4 )) 

i=l j=3 

xe^^-^e-^'^-^rK^-kKKs'h-k'), (14) 
where = (si + s 2 )/2 and s' h = (s 3 + S4)/2. Repeatedly using the distant observer approximation, 

J d Sl J ds 3 V(s 1 ,k)^(s 3 ,k')(P(s 1 )F(s 3 ))e lk - sl - 4k '- S3 

~ J dsn(s)V(s,k)V(s,k') e i;s - (k - k,) [P(^~, + ^] ■ (15) 

we finally have 



(AV l (k)AVl(k'))~AP?(k)S {3) (k~k , ), (16) 
where 

_ 0.0 , 3 Jrfsn(s)V(s,k) 4 [P(k, |sj) + (1 + «)/n(s)] 2 [£,(k • s)] 2 
A ^ W ~ 2(2?F) [/ds'n(s') 2 V(s',k) 2 ]2 ■ ^) 

We define the estimator of the multipole moments of the anisotropic power spectrum by 

Pi(*)=aW rfk ?M k )> (18) 

where AT4 denotes the volume of the shell in the Fourier space. The ensemble average of Vi(k) is given by 

IVtlrW 1 f ^ J^)V(s ; k) 2 P(k,|s|)£ ; (k-s) 
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Note that Vi(k) corresponds to the multipole coefficient Pi(k,z) in equation (1). Vi(k) is generalized so as to incorporate 
the redshift evolution. In the limit that the sources are confined to a very narrow range of redshifts z, (Vi(k)) reduces 
to Pi(k,z). 

The variance of Vi(k) is obtained by evaluating 
(A^(fc) 2 )^([^(fc)-(^(fc)>] 2 > 

= aW dk / rfk'(A7Mk)A7Mk')>, (20) 

^ V k JAV k JAV k 



which reduces to 

(AV l (k) 2 } = 2^Q?(s,fc), (21) 
where we have defined 



1 f , Jrfsn(s)V(s,k) 4 [P(k, |s|) + l/n(s)]'[£,(k ■ s)] 2 
Q < (S ' fc) " AV k J AVk dk [/dB'n(B')W,k) 2 ] 2 
Here we have assumed a«l. 
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Now we find that the following weight factor gives the stationary solution 
^ k)= l + ^)P(k,|s|) (23) 



for the variation 
S(AVi(k) 2 ) 



0. (24) 



<ty(s,k) 

In this case find the minimum error 



AV k 2j_, P /ds'fi(s') 2 [l + n(s')P(fc,M,|s'l)] 

where we used the notation fj, = k- s, ^(Sjk) = ^(s, fc,/x) and P(k, |s|) = P(k,fj,, |s|). In the case nP(k, |s|) <C 1, we may 
set ^(s,k) = 1 and we have 

(2, + m p l{k f) - <A W )*> - S^J^, (26) 
where we assumed the Legendre polynomials are normalized by 

dfiC{fjL)£ v (fjL) = ^—6 w . (27) 



Thus the error of the Z-th moment is smaller than in the monopole spectrum by y2l + 1. 

Here we summarize the difference between our method and that developed by Cole et al. (1994). Using our notation, 
essentially, they first estimate 

/ dgi / ds 2 ^( Sl , k)^(s 2 , k)F( Sl )F(s 2 )e- k -^-^) _ 

Jrfsn 2 (s)^( S ,k)2 (2 * } 

which is same as the right hand side of equation (7) but without Ci(sh -k). After averaging i?(k) over the angle ip 
(Figure 1), and subtracting the shotnoise contribution, the anisotropic power spectrum P(k,[i) can be obtained. Then 
the quadrupolc component can be obtained from P(k,[i) using Ci(jx). In oder to define the direction of the line of 
sight, in their method, a Gaussian window function is introduced, and the above procedure is repeated for the sample 
with different windows and the result is averaged. Essentially, in their method, the sample is divided into subsamplcs 
in order to clearly introduce the quantity /i associated with a line of sight direction. Thus the method by Cole 
et al.(1994) first requires computation of the anisotropic power spectrum, whilst our method does not. The estimation 
of the multipole spectrum can be performed completely parallel with that for the estimation of the monopole spectrum. 
However, in principle, both the methods are equivalent in the limit of the distant observer approximation (excepting 
the weight factor). 

2 Compare with the result in (Yamamoto 2003), in which the weighting factor was determined under the constraint that it has no angular 
dependence -0(s,k) = i/j(s, fe). 
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Next let us consider the estimator for the discrete density field of the object catalog. In this case we replace £j(sfc-k) 
in equation (7) with Ci(s\ ■ k), which is valid as long as the distant observer approximation is applicable, and then we 
may write 



N N s 



R t (k) = A- 1 \J2 1>{*ii , k)e i;k ' Sl1 & (s 4l • k) - aJ2^n , k)e* k ^ & (s n ■ k) 




^V(s l2 ,k)e 4k - S - - a ^V(s, 2 ,kK k - s - , (29) 



.12 



where N (N s ) is the number of objects of the real (synthetic) catalog, respectively, and we have defined A = 
J dsn(s) 2 ip(s,]c) 2 . In a similar way, the shotnoise term may be written as 

TV 

Si(k) = A-^l + a)^^^) 2 ^ ■ k), (30) 

i 

using the real catalog, or 

5j(k) = A-^l + a)a^(s,,k) 2 £ i (s, • k), (31) 

3 

using the synthetic catalog. The other possible expression is 

(N N s \ 

2V(si,k) 2 A(si-k) + a 2 2^( s i' k ) a/: '( S J- k ) ' ( 32 ) 

using both the real and synthetic catalogs, as used in Cole et al. (2005). The choice of the shotnoise estimation does 
not alter our results in the present paper. 

3. Application to the 2QZ survey 

In this section, we apply the method to the 2dF QSO sample. In general, the constraint on the cosmological 
parameters from the QSO sample is not very tight, however, it is useful to obtain unique information about the high 
rcdshift universe. The clustering of the sample has been investigated by the 2QZ group extensively (e.g., Outram 
ct al. 2003, Croom et al. 2005 and references therein). The result from the clustering analysis is quite consistent with 
the prediction of the concordance A cold dark matter (CDM) model built up from consideration of other cosmological 
data (e.g., Spcrgcl ct al. 2003, Cole et al. 2005, Tegmark ct al. 2004, Riess et al. 1998, Perlmuttcr et al. 1999). The 
redshift-space distortion in the clustering of the 2QZ sample has been investigated by Outram et al. (2004). In their 
work, the anisotropic power spectrum P(k^,k ± ) is measured, and the value of (3(= fl^/b) is determined, where fcy and 
k ± are the wave numbers parallel and perpendicular to the line of sight direction, respectively. In contrast to their 
approach we focus on the quadrupole spectrum in the clustering of the 2QZ sample. 

The present work is an extension of previous work investigating the monopole power spectrum with the 2QZ sample 
(Yamamoto 2004). Here we briefly review the survey details for self-containment. The 2QZ survey covers two areas 
of 5 x 75 deg 2 , one in the south Galactic cap (SGC) and the other in the north Galactic cap (NGC) in the redshift 
range of less that 3. The survey area is defined by the equatorial coordinates from a = 21 h 40 to a = 3 h 15 and 
-32.5° < 6 < -27.5° in the SGC, and 9 h 50 < a < 14 h 50 and -2.5° < 6 < 2.5° in the NGC, respectively. 

We use the final catalog reported by the 2QZ survey group, which is available through the home page 
http://www.2dfquasar.org/. Details of the QSO selection of the catalog is also described in the reference Croom 
ct al. 2004, in which the photometric and spectroscopic incompleteness is discussed, as well as the angular selection 
function. The coverage incompleteness is given by the hole information which is publicly available through the home 
page. We generate the random sample through a random process mapping angular position of a QSO to different 
position with fixing the redshift with an equal probability on each survey area. Then we use 10713 and 8442 QSOs in 
the SGC and the NGC, respectively, in the range of redshift 0.2 < z < 2.2 incorporating the hole information. Outram 
et al.(2004) used the QSOs in the range of redshift 0.3 < z < 2.2, which is slightly different from our choice. This 
difference does not alter our result because the number of the QSOs in the range 0.2 < z < 0.3 is less than one percent 
of the total number of the QSOs used here. 

Next we explain our theoretical modeling of the power spectrum. In a redshift survey, the rcdshift z is the indicator 
of the distance. Therefore we need to assume a distance-redshift relation s = |s| = s(z) to plot a map of objects. The 
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power spectrum depends on this choice of the radial coordinate of the map s 
For the A CDM model the comoving distance is given by 

f z dz' 



ff ov /fi m (i+2') 3 + i-a™ 

100/ikm/s/Mpc is the Hubble parameter. 



s(z) due to the geometric distortion. 



(33) 



For our fiducial model we adopt the fiat ACDM model with 



where Ho 

= 0.3 Thus our fiducial model is s(z) = r(z, 0.3). 
Our theoretical model properly incorporates the light-cone effect (redshift-cvolution effect of the clustering), the 
geometric distortion, the linear distortion, the Finger of God effects, including the nonlinear evolution of the density 
perturbations. We model the QSO power spectrum as (e.g., Suto et al. 2000b, Yamamoto 2002; 2003) 



P(k,s(z))=P(k, t i,s(z)) 



s(z) 2 ds(z) 
r(z) 2 dr(z) 



Pqso q 



ds(z) 
k/J, , , r ,q± 



dr(z) 



k^l 



■£(£) 
r(z) 



with 



PQSo(qn,q±)=b(zy[l + 



d\nDi(z)/d\na(z) q { \ 



b(z) 



for the nonlinear modeling, and 

dlnD 1 (z)/d\iaa(z) q 2 



PQSo{q\\,q±) =b{z) 



1 



b(z 

2 



q 



P 



p L 



Nonli 



\q,z)D[q^ P {z)]V{5z) 



\q,z)V{8z) 



(34) 



(35) 



(36) 



for the linear modeling, where q- — q {l +q ± , D\(z) is the linear growth factor normalized to D±(z = 0) = 1, ^mass^fe- 2 ) 
( • p mas n s lin ° ar (9^)) & the linear (nonlinear) mass power spectrum, D[q ll ap(z)] is the damping factor due to the Finger of 
God effect, and T>(Sz) represents the damping factor due to the error in rcdshift measurement (see equation (37)). We 
model the exponential distribution function for the pairwise peculiar velocity, using an approximate formula for the 
mean square velocity dispersion at the large separation determined through the cosmic energy equation (Mo et al. 1997, 
Magira et al. 2000, Suto et al. 2000a). In particular for the QSO sample, the error in measuring the redshift is not 
negligible, which causes an apparent velocity dispersion and additional redshift-space distortion. We incorporate this 
effect in modeling the power spectrum by multiplying the spectrum with the damping factor 



V(5z) = cxp 



1.2 



<5z 2 > 



(37) 



where fcy is the comoving wave number parallel to the line of sight direction, 8z is the variance in the redshift 
measurement error for which we adopt Sz = 0.0014(1 + z) (Croom et al. 2005). 

In the modeling of the bias, we consider the scale independent bias model of Fry (1996) 



b(z) 



1 



h-1 

D 1 (zY 



(38) 



where bo is the constant bias parameter. 

Figure 2 plots our theoretical curves for Po(k) (left panel) and Viik) /Vo(k) (right panel), to show which effect is 
important in the power spectrum. The upper (lower) panel assumes bo = 1 (bo = 2) in the bias model (38). Here we 
adopt the other cosmological parameters Sl m = 0.28, fib = 0.045, h = 0.7, erg = 0.9 and n = 1 motivated the WMAP 
result (Spergel et al. 2003). In this figure the dashed curve is the linear modeling (36), while the solid curve is the 
nonlinear modeling (35). Concerning the monopole spectrum Vo> the linear model is very similar to the nonlinear 
model. This occurs by a cancellation effect of the two nonlinear effects; One is the increase of the amplitude of the 
real-space power spectrum at large k. The other is the decrease of the amplitude due to the damping factor of the 
Finger of God effect D[q ll ap(z)] and the error in redshift measurement T>(dz). Concerning the quadrupole divided by 
the monopole, V2/V0, the behavior at small k comes from the linear distortion effect, while the behavior at large k 
comes from the damping factor T>(Sz) (dashed curve). In the absence of V(Sz), V2/V0 is constant for k in the linear 
modeling. In the nonlinear model of V2/V0 (solid curve), the further decrease appears due to the Finger of God effect, 
thought the effect is not very significant. This suggests that the error in redshift measurement is the significant effect 
in the present sample for V2 ■ The effect of the geometric distortion is negligibly small in our model because the fiducial 
model to plot map and the theoretical model is very close. Namely, s(z) and r(z) are not very different. 

Then we followed the prescription explained in the previous section for the multipole moments of the power spectrum. 
We generated the random catalog with the parameter choice a = 1/15. Figure 3 plots the monopole power spectrum 
(open squares on the upper panels) and the quadrupole spectrum divided by the monopole spectrum (open squares on 
the lower panels). The left and right panels show the results with the SGC and NGC samples, respectively. The error 
bars are estimated from the expression (26). The solid curve in Figure 3 is the theoretical curve of the ACDM model 
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with the cosmological parameters same as those of Figure 2, and bo = 1.5 the best fit value of the bias parameter. For 
small wave numbers less than k < 0.03/iMpc -1 , the measured quadrupole has negative sign. 

Here let us briefly compare our approach with that of Outram ct al. (2004). As mentioned before, they measured 
the anisotropic power spectrum P(k ll ,k ± ), while we have measured the quadrupole moment l D 2{k). In their analysis, 
in order to estimate P(k n ,k ± ), each 5 x 75 deg 2 strip was divided into 8 regions. And for each of these small region, 
fen and k ± were defined by applying the distant observer approximation. Thus P(k ]]7 k ± ) were estimated separately 
in 16 individual 5 x 10 deg 2 regions. We have also estimated the quadrupole spectrum in these divided samples, but 
the method clearly failed to detect the quadrupole spectrum. As the area of the individual region becomes small, the 
estimated values of the quadrupole spectrum become negative. This suggests that the quadrupole spectrum with our 
method is sensitive to the shape and the boundary of the survey area of the sample. 

Earlier work investigated the constraints from 2QZ on the cosmological parameters f2 m , and the equation of state 
of the dark energy vo, by fitting to the monopolc power spectrum (Yamamoto 2002, 2004). In the present work, by 
adopting the above theoretical model, we compare the measured quadrupole spectrum with the theoretical prediction. 
Using the bias model (38), for this comparison we define 



where T'i(k i )° h is the value of the observed power spectrum at k i} A7';(fc i ) is the variance of errors in Figure 3, and 
Pi(ki) th is the theoretical spectrum predictions. Figure 4 shows % 2 a s a function of the bias parameter bo- The left 
and right panels show the results for the SGC and NGC samples, respectively. Here we fixed the other cosmological 
parameters to be the same as those in Figure 2. The degree of freedom is 12. In each panel, the dashed curve is the 
result from To, the best fit value of bo is 1.4 (1.5) for the SGC (NGC) sample. On the other hand, the solid curve is 
the result from the quadrupole, Vn- Because the error of Vi is larger than from Vq, the constraint on bo from V2 is 
weak. However, we see the minimum of \ 2 of the solid curve is located around the minimum of the dashed curve in 
the left panel. These are consistent with the best fit value of bo — 1.64 found by the 2QZ group though the assumed 
cosmological model is slightly different from ours (Croom et al. 2005). 

For the right panel (NGC), the agreement is not good compared with the SGC sample, but we see roughly the 
same behavior, although there are some troubling discrepancies that deserve further study. We may conclude that the 
quadrupole signal is detected in our analysis, though the error is rather large. As the reason of the poor detection, 
we infer a contamination of the systematic error due to the survey geometry, i.e., the shape and the boundary effect 
of the inhomogencous survey region. 

Here we discuss the validity of bias model. As is mentioned by Croom ct al.(2001;2005), the redshift-evolution of 
bias (38) does not fit observational data very much. They have given the best fit formula (Croom et al.2005). 



The dashed curve in Figure 3 plots the theoretical one using this fit formula of the bias evolution. We see the good 
agreement between the dashed curve and the solid curve. Note that the solid curve is the theoretical curve using the 
bias model (38). This agreement comes from the fact that our theoretical spectrum depends only on a mean value of 
the bias. 

Finally in this section, we briefly mention the validity of the distant observer approximation. Cole et al.(1994) 
discussed various effects which affect the proper /3-factor estimation from the redshift-space distortion. A possible 
relevant effect here is the large opening angle on the survey. Cole et al.(1994) found that a correction might be needed 
to be applied for angles greater than 30 degree. However, the range of the redshift of the 2QZ sample is large. At a 
mean redshift z = 1.2, the wave number k = 2n/X is less than 0.005 hMpc^ 1 for A = s(z)sm(8) with 8 = 30 degree. Thus 
the distant observer approximation is good in our situation. On the other hand for 6 = 5 degree, k ~ 0.03 /iMpc -1 , 
therefore deviations in the observed quadrupole spectrum away from the theoretical prediction may appear around 
these wavenumbers. 

4. Conclusions 

In summary we have investigated the quadrupole of the redshift-space power spectrum, P2(k), of the 2dF QSO 
survey. First, we have developed an algorithm which estimates the quadrupole and monopole spectra in a unified 
manner. The error in a measurement of the quadrupole spectrum is given analytically. We have detected the signal 
of the redshift distortions in the quadrupole spectrum, though the errors are still substantial. 

When we used the sample divided into several subsamples, the agreement of the estimated spectrum with the 
theoretical model became worse. However, the resulting power spectrum is consistent with the theoretical prediction 
and with previous analysis of the 2QZ survey. Concerning the bias, however, as discussed by Croom et al. (2001;2005), 
the Fry bias model is not the best fitting model for the evolution of the QSO bias. To check a more realistic model of 




(39) 



b(z) = 0.53 + 0.289(1 + z) 2 . 



(40) 
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the bias, we used the formula (40) in the theoretical modeling, which best fits the data from Croom et al.(2005). We 
found that this theoretical curve is almost the same as that of the Fry's bias model. This is because our theoretical 
spectrum depends only on the mean value of the bias. 

Constraining the nature of dark energy is one of the ultimate purposes of such a clustering analysis. However, as 
demonstrated in (Yamamoto 2004), the 2dF QSO sample cannot provide a stringent constraint on the equation of 
state of the dark energy. The reason for this is the large shotnoise in the QSO distribution since the mean number 
density is too small. The application of the algorithm tested here to other, larger, galaxy samples such as SDSS will 
be very interesting (see e.g., Eisenstein et al. 2005, Cole et al. 2005, Yahata et al. 2005) and is left to future work. 
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Fig. 1. A Sketch for the definition of the variables. 
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Fig. 2. Theoretical prediction for the monopole power spectrum Vo{k) (left panel) and quadrupole V2(k) /'Po(fc) (right panel). The 
dashed curve is the linear modeling (36), while the solid curve is the nonlinear modeling (35). The upper (lower) panels assume 
bo = I (bo = 2) in the bias formula (38). Theoretical model is the ACDM model with the cosmological parameters H m = 0.28, 
Q b = 0.045, h = 0.7, ct 8 = 0.9 and n = 1. 




Fig. 3. The monopole power spectrum Vo{k) (upper panels) and quadrupole divided by monopole V2(k)/Vo{k) (lower panels) from 
the 2dF QSO sample. The left (right) panels are from the SGC (NGC) sample, respectively. The curves assume the ACDM model 
with the cosmological parameters Q m = 0.28, Qt, = 0.045, h = 0.7, erg = 0.9 and n = 1. Concerning the bias model, the solid curve 
assumes (38) with 6q = 1-5, but the dashed curve uses (40). 
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Fig. 4. The \ 2 f° r T'o(fc) (dashed curve) and "Ps^k) (solid curve) as a function of the bias parameter bo. The left (right) panels are 
for the South Galactic Cap (North Galactic Cap) sample, respectively. 



